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Abstract 

The angular correlation of the electrons emitted in the neutrinoless double beta decay {Ov20) is 
presented using a general Lorentz invariant effective Lagrangian for the leptonic and hadronic charged 
weak currents. We show that the coefficient K in the angular correlation dF/dcosO oc (f — KcosO) is 
essentially independent of the nuclear matrix element models and present its numerical values for the five 
nuclei of interest C^Ge, *^Se, ^°°Mo, -^^°Te, and ^^''Xe), assuming that the Oi'2/3-decays in these nuclei 
are induced solely by a light Majorana neutrino, vm. This coefficient varies between K = 0.81 (for the 
''^Ge nucleus) and K = 0.88 (for the *^Se and '^''"Mo nuclei), calculated taking into account the effects 
from the nucleon recoil, the S and P-waves for the outgoing electrons and the electron mass. Deviation 
of K from its values derived here would indicate the presence of New Physics (NP) in addition to a light 
Majorana neutrino, and we work out the angular coefficients in several um + NP scenarios for the '^''Ge 
nucleus. As an illustration of the correlations among the 01/2/3 observables (half-life T1/2, the coefficient 
K, and the effective Majorana neutrino mass |(?n)|) and the parameters of the underlying NP model, we 
analyze the left-right symmetric models, taking into account current phenomenological bounds on the 
right-handed Wjj-boson mass and the left-right mixing parameter (. 
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1 Introduction 



It is now established beyond any doubt that the observed neutrinos have tiny but non-zero masses and they 
mix with each other, with both of these features foUowing from the observation of the atmospheric and 
solar neutrino oscillations and from the long baseline neutrino oscillation experiments [T]. Theoretically, it 
is largely anticipated that the neutrinos are Majorana particles. Experimental evidence for the neutrinoless 
double beta decay (0i^2/?) would deliver a conclusive confirmation of the Majorana nature of neutrinos, 
establishing the existence of physics beyond the standard model. This is the overriding interest in carrying 
out these experiments and in the related phenomenology [2] . 

We recall that Oz^2/3-decays are forbidden in the standard model (SM) by lepton number (LN) conser- 
vation, which is a consequence of the renormalizability of the SM. However, being the low energy limit of 
a more general theory, an extended version of the SM could contain nonrenormalizable terms (tiny to be 
compatible with experiments), in particular, terms that violate LN and allow the 0v2P decay. Probable 
mechanisms of LN violation may include exchanges by: Majorana neutrinos vms [3j [H [5] (the preferred 
mechanism after the observation of neutrino oscillations [T]), SUSY particles [51 [71 [51 [HI UHl EI] , scalar bilin- 
ears (SBs) [12], e.g. doubly charged dileptons (the component <^^^ of the SU{2)l triplet Higgs scalar etc.), 
leptoquarks (LQs) [12|, right-handed Wr bosons [5l [Tlj etc. From these particles light i/s are much lighter 
than the electron and others are much heavier than the proton. Therefore, there are two possible classes of 
mechanisms for the 0v2/3 decay. With the light z/s in the intermediate state the mechanism is called long 
range and otherwise it is referred to as the short range mechanism. For both these classes, the separation 
of the lepton physics from the hadron physics takes place [15], which simplifies calculations. According to 
the Schcchtcr-Valle theorem [16], any mechanism inducing the 0i^2/? decay produces an effective Majorana 
mass for the neutrino, which must therefore contribute to this decay. These various contributions will have 
to be disentangled to extract information from the Qi^2f3 decay on the characteristics of the sources of LN 
violation, in particular, on the neutrino masses and mixing. Measurements of the neutrinoless double beta 
decay in different nuclei will help in determining the underlying physics mechanism [171 118j . 

Our aim in this paper is to examine the possibility to discriminate among the various possible mechanisms 
contributing to the Oz^2/3-decays using the information on the angular correlation of the final electrons in 
the process Ni{A, Z) Nf{A, Z + 2) + e~ + e~ . A preliminary study along these lines was published by 
us in 2006 [19], with admittedly simplified treatment neglecting the nucleon recoil and the P-wave effects 
in the outgoing electron wave function. We rectify these shortcomings and provide in this paper a detailed 
account of the improved treatment. Restricting ourselves to the long-range mechanism, treating the electrons 
relativistically but with non-relativistic nucleons, we derive the angular correlation between the electrons 
using the general Lorentz invariant effective Lagrangian involving the leptonic and hadronic charged weak 
currents. Generally, this angular correlation can be expressed as dV/d cos 6 ~ 1 — Kcos9, where 6 is the 
angle between the electron momenta in the rest frame of the parent nucleus. Expressing K — B/A, with 
— 1 < JsT < 1, we derive the analytic expressions for A and B for the effective Lagrangian characterized by 
the coefficients e^j encoding the standard, {V — A) {V — A), and new physics contributions (see Eq. (1)). 
Essential steps of these derivations are presented in section 2. The analytic expressions derived here confirm 
the earlier detailed derivations by Doi et al. [S] , and we specify where the treatment presented here transcends 
the earlier work. Specific cases are relegated to Appendix A (for the decays involving scalar nonstandard 
terms). Appendix B (for the vector nonstandard terms), and Appendix C (for the tensor nonstandard terms). 
We hope to return to the discussion of including the short-range mechanism, neglected in this paper, in future 
work. 

Numerical analysis of the electron angular correlation is presented in section 3, and the coefficient K for 
the various underlying mechanisms in Oz^2/3-decays are worked out. In particular, numerical values of K for 
the five nuclei of current experimental interest: ^^Ge, ^^Se, ^"°Mo, ^'^"Te, and ^'^^Xe are presented for the 
light Majorana neutrino I'm case. Their values range from K — 0.81 (for the ^^Ge nucleus) and K — 0.88 
(for the ^^Se and ^*"^Mo nuclei). To study the uncertainty in the nuclear matrix elements, we have employed 
the so-called QRPA model with and without the p-n pairing for the ^^Gc nucleus [ID], and a more modern 
QRPA model, fixing the particle-particle pairing strength [21]. While the uncertainty due to the nuclear 
matrix element model is quite marked for 2^/2 in some cases, we show that it is rather modest for K, not 
exceeding 10% for the models discussed here. For the + NP scenarios, we remark that the nonstandard 
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coefRcients e^ipyi' ^Th' ^"^"^ ^t^ '^'^ '^'-'^ change the value of the angular coefficient K. The contribution of the 
scalar nonstandard term from the ef coefficients is found to be numerically small. So, what concerns the 
angular correlation, we have essentially three distinct scenarios: (i) Standard {vm), (ii) R-parity violating 
SUSY [vm + ^Tfl)' ^"^^ (iii) left-right-symmetric models {vm + ej/!]!^)- Numerical analysis of the coefficient 
K in the extended :/m + NP scenario is carried out for the decay of the ''^Ge nucleus using the nuclear matrix 
element model already specified. 

We take a closer look at the underlying physics behind the coefficients e^^^ in section 4. These coefficients 
appear in the context of the left-right symmetric models which are theoretically well motivated |22| . Also, 
the corresponding nuclear matrix elements are available in the literature. Making use of them, we work out 
the correlations among the angular coefficient K , the half- life T1/2 and either the mass of the right-handed 
Wr boson, tuwr, or the W boson's mixing angle C, taking into account the current bounds on the various 
parameters. Results are presented in Figs. 1-4. The differential distribution dT /d cob for the Qv2j3 decay 
of the ^^Ge nucleus is shown in Fig. 5 for some representative values of |(m)| for tuwr = 1, 1-5 TeV and 
for an infinitely heavy mwn ■ It is seen that the effect of the right-handed W^-boson is more marked in the 
angular correlation for smaller values of | (to) | . 

2 Angular correlation for the long range mechanism of decay 
2.1 General effective Lagrangian 

For the decay mediated by light i^A/s, the most general effective Lagrangian is the Lorentz invariant combi- 
nation of the leptonic ja and the hadronic Jq currents of definite tensor structure and chirality [231 124] 

C = ^^[(C/e. + el-_Xdj^'_^4_A., + Y^'^^^yi + H-c] , (1) 

^ «,/3 

where the hadronic and leptonic currents are defined as: = uOad and = eOpVi\ the leptonic currents 
contain neutrino mass eigenstates and the index i runs over the light eigenstates. Here and thereafter, a 
summation over the repeated indices is assumed; a,f3—V^A,S^P,TL^ii^ [Ot^ = 2a'^Tp, a'^" = | [7^,7"], 
Pp = (l=P75)/2 is the projector, p ~ L, R); the prime indicates the summation over all the Lorentz invariant 
contributions, except ior a ^ /3 = V — A, U^i is the PMNS mixing matrix 25J and Vud is the CKM 
matrix element [1 . Note that in Eq. ^ the currents have been scaled relative to the strength of the usual 
V — A interaction with Gp being the Fermi coupling constant. The coefficients e^^ encode new physics, 
parametrizing deviations of the Lagrangian from the standard V — A current-current form and mixing of the 
non-SM neutrinos. 

In discussing the extension of the SM for the 0v2(3 decay, Ref. [S] considered explicitly only nonstandard 
terms with 

9v 9v 
Implicitly, also the contributions encoded by the coefficients e^I^ j are discussed arising from the non-SM 
contribution to Uei in SU{2)l x SU{2)r x [/(I) models with mirror leptons (see Ref. [5], Eq. (A. 2. 17)). Here 
V, U' and V are the 3x3 blocks of mixing matrices for non-SM neutrinos, e.g., for the usual SU{2)l x 
SU{2)fi X C/(l) model V describes the lepton mixing for neutrinos from right-handed lepton doublets; for 
SU{2)l X SU{2)ii X U{1) model with mirror leptons [26] U' (V) describes the lepton mixing for mirror 
left(right)-handed neutrinos [5] etc. The form factors gv and g'y are expressed through the mixing angles 
for left- and right-handed quarks. Thus, gv = cos 9c = Vud and g'y = e*^ cos^^, with 9c being the Cabibbo 
angle, 9'(-, is its right-handed mixing analogoue, and the CP violating phase S arises in these models due to 
both the mixing of right-handed quarks and the mixing of left- and right-handed gauge bosons (see Ref. [5] , 
Eq. (3.1.11)). The parameters k, 77, and A characterize the strength of nonstandard effects. Below, we give 
some illustrative examples relating the coefficients eyl^ j, ^v±A i ^^'^ the particle masses, couplings and the 
mixing parameters in the underlying theoretical models. 

In the R-parity- violating (RPV) SUSY accompanying the neutrino exchange mechanism [6l[7l[8l[9lfT0|,lllj. 
SUSY particles (sleptons, squarks) are present in one of the two effective 4-fermion vertices. (The other vertex 
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contains the usual Wl boson.) The nonzero parameters are 

1 . 



-V-A,: 
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where the index n runs over e, ^, r (1, 2, 3), and the RPV Minimal Supersymmetric Model (MSSM) 
parameters r/s depend on the couplings of the RPV MSSM superpotential, the masses of the squarks and 
the sleptons, the mixings among the squarks and among the sleptons. Concentrating on the dominant 
contributions ef^pj and e^^ ■ (as the others are helicity-suppressed) , one can express V^q-^LP. ^^^"^ '^(ijLfl 
follows [TD] 
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where k is the generation index, 0^^^ and 0^^.^ are the squark and slepton mixing angles, respectively, mj^ 
and nif are the sfermion mass eigenvalues, and Xijk and A^ j. are the RPV-couplings in the superpotential. 



For the mechanism with LQs in one of the effective vertices [T^ , the nonzero coefhcients are 
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where 



(5) 



(6) 



the parameters es{v)i Q^s^y)' ^^'^ '^^S{v) depend on the couplings of the renormalizable LQ-quark-lepton 
interactions consistent with the SM gauge symmetry, the mixing parameters and the common mass scale 
Ms(v) of the scalar (vector) LQs [77] . 

The nonzero for the discussed models are collected in Table 1. 

Table 1: Nonzero coefficients for various models. 



AR) 



Model 


Nonzero es 


with Wrs 


y-A y+A 

^V+A' ''V+A 


RPV SUSY 


^s+p^ y^A^ Tr 

'^S+P^ ^V-A' ^Tr 


with LQs 


s+p y+A 

^StP' ''V+A 



The upper bounds on some of the parameters ([6|) from the Heidelberg-Moscow experiment were 
derived in Ref. [28j using the S'-wave approximation for the electrons, considering nucleon recoil terms and 
only one nonzero parameter e^^ in the Lagrangian ([1]) at a time. 



The coefficients e^^ entering the Lagrangian ^ can be expressed as 



(7) 



where U^f'^'^ are mixing parameters for non-SM neutrinos (see, e.g., Eq. As this Lagrangian describes 

also ordinary /3-decays (without LN violation), the coefficients are constrained by the existing data on 
precision measurements in allowed nuclear beta decays, including neutron decay [29] . For example, from 
these data we obtain the conservative bound 



|e^+^| <7x 10- 



(8) 
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From Eqs. ©, (O and the bound |e 

nonstandard mixing is small: 



V+A\ 

v+a\ 



< 7.9 X 10 ^ (see section 3.2) we can assume that the 



\Ue^Ve^\ < 10~^ V,, = U. 



(V+Ay+A) 



(9) 



2.2 Methods and approximations 

We have calculated the leading order in the Fermi constant taking into account the leading contribution of 
the parameters to the decay matrix elements using the approximation of the relativistic electrons and 
non-relativistic nucleons. The wavefunction of an electron with the asymptotic momentum p and the spin 
projection s can be expanded in terms of spherical waves as O I30j 



eps(r) 



We take into account the iS'1/2 and the P1/2 waves for the outgoing electrons: 



Pl/2 



(r) 



J-lXs 

ho- ■ VXs 



(10) 

(11) 

(12) 



with r = r/r, p = p/p and the two component spinor Xs- We use the approximate radial wave functions [5] 
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f-i 



r 11 
= ±A^iC±(e)-^, = -aZ + -{e ± me)R, 



(13) 
(14) 
(15) 



including the finite de Broglie wave length correction (FBWC) for the 5*1/2 wave. Here R is the nuclear 
radius, e is the electron energy and a is the fine structure constant. For the normalization constants A±i 
we use the approximate Eq. ((45|l (see below). 

The nucleon matrix elements of the color singlet quark currents are [U [SU [32l [33] 



{P{k')\u{lTl,)d\N{k))^{l^{k') \Ff{q^)TFP{q^)^,] T+m, 



{Pik')\ua'^-'ilTl5)d\Nik)) =i^ik') 



9v{q^h^ T 3^(9^)7^75 - mi{q^)^ — - ± gpiq^hsq^ 



2m,- 



T+'ipik), 



(16) 

r+V(fc),(17) 
(18) 



where 



-^(3) ^(3) 
jM- ^ Tf '(g2)aP- + -^ij^^q" - -f'^qn + -^{cj'^''qpq'' ' o'^Pqpq'^), (19) 

(20) 



is a nucleon isodoublet. 

The non-relativistic structure of the nucleon currents in the impulse approximation is derived using Refs 
j32|l34]. see Appendices A, B, and C. We have calculated the nucleon recoil terms including the recoil terms 
due to the pseudoscalar form factor. 
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Table 2: Expressions for A in Eqs. ([^T|l and for the stated choice of e^. 
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Table 3: Expressions for B in Eq. ([221) for the stated choice of e[ 
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^V-A 


Bo + AD,M^,^Zi\co2+^D^\fi)^Z^\^ 
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^V+A 
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Bo 
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2.3 Electron angular correlation 

Taking into account the dominant terms introduced in the Appendices A, B, and C in the closure approxi- 
mation [5j we obtain the differential width in cos0 for the 0^{A, Z) — >0+(^, Z + 2)e^e^ transitions: 

dr In 2 , , n , , , , , 

-\Mgt\ A{1- K cos 9), (21) 



d cos 9 2 

where 9 is the angle between the electron momenta in the rest frame of the parent nucleus and the angular 
correlation coefficient is 

X = ^ , -1< X < 1. (22) 

The Gamow-Teller nuclear matrix element Mqt is defined in Eq. (|5ip below. 

The expressions for A and B for different choices of e£, with only one coefficient considered at a time, 
are shown in Tables 2 and 3. 



In these tables 



and 



i = cos^Z-i, Si = sinV'i (23) 
{■m)/me, ijlI^ ml/me, (24) 
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with the standard effective Majorana mass (m) = U'^i'mi and the nonstandard ones: 

i i i 

The quantities A and B for all zero are 

^o = Ci|m|', Bo = Di|/x|2 

and the relative phases are 

^01 = arg((/x)/x|^;^*), Vo2 = arg((/x)/xj:^l^*), 

^1 = arg((M)e^+t), ^2 - arg((M)e|;+t), 

-003 = f^Tg{{fi)fil.-^*), Vo4 = arg((/x)/x|l^*), 

i)3 = arg((/z)e|+^*), 1/^4 = arg((/z)e|+^*), 

V'06 = arg((/i)/z^^*), 

1^5 = arg((/i)e^^*), Ve = arg((/i)e^«*). 

The coefficients Ci and Q ' in Table 2 are 

Co = {xl-l)Aoi, 
Ci = ixF-l?Aou 

C2 = {XF - l)(X2-^03 - Xl+^04), 

Cs = -(XF-l)(X2+^03-Xl-^04-Xp^05 + xk^06), 

2 1 

Ci = X2-A02 - -Xl + X2-^03 + gXl+^04, 

2 1 

C5 = X2+^02 - gXl-X2+-4o3 + gXl-^04 + Xp-4o8 - x'pXk^O? +Xfl-4o9 
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The coefficients Di and d'-^^''^^ entering in Table 3 are: 
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where the integrated phase space factors are 
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with the phase space element dQo„ defined as follows: 

dQov = TO^^|pi||p2|eie2^(ei +S2 + Ef - Ei)deids2d{pi ■ P2) . 

(SPT) i S P T) 

The constant aoi/ and the kinematic factors oofc, ' , 6ofc and 6qj.' ' ' entering above are defined 
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where £21 = £2 — £1 is the difference in the electron energy. The characteristic features of the Pi/ 
expressed as 

C = SaZ+{ei +e2)R 

and the Coulomb corrections appear as the following combinations 

a± = |q;_i_iP ± |aiip, /3± = |ai-ip ± 
7+ = 2Re(aiiQ;li_i), 7_ = 2lu\{aiia*_i_i) , 
(5_l_ — 2'Re{a-iia\_i) ^ (5_ = 2Im(Q;_iia*_]^), 

with = Ai{£2)Aj{ei). 

For the normalization constants in the approximation including terms up to {aZY' [5] 
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we have 
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Note that using Eq. (|T7)l the expressions for B from Table 3 are reduced to the form shown in Table 4. 
Table 4: Expressions for B in Eq. for the stated choice of for the A±i from Eq. 
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In the definitions of d and Di we use some combinations of nuclear parameters similar to the ones in 
Ref. 0. Thus, 

X2± = XGTuj ± XFu - ^XiTi Xi± = (x'gt - 6Xt) ± 3xf; 



, 9v 

XF = 



Xfe 



Mf . 

Mgt' 
Mfe 



Xfc 



gv Mk 



M, 



gA Mgt 
k^T, GT, RCa, RTy, 



k:=^ R, RT: 



Xf 



(3) 



— x-f; Xfo 



GT 



(3) 



g_v\ 



M 



FO 



Xpo — ,^ ; 

gA gA Mgt 



F 



gv 

(3) 



gv V5A/ Mgt' 



Xl = — Xfe, k = R, RT, RG^, RT,; xl 

gA 



) Ml , 
— k 

gA Mgt 



(49) 



where the index F refers to Fermi, GT to Gamow-Teller, T to tensor, P to the P-wave effect and R to 
the recoil effect. If x has prime or the index uj than the same has the according matrix element in the 
numerator. The nuclear matrix elements defined below contain the operator = (ti + iT2)"/2 converting 
the a-th neutron into the a-th proton, and the initial (final) nuclear state are denoted by |0^) ((Oj |) 



MF = Y.^O+\\Y,h+{rab.EN)TlrX\\Qt). 

N a^b 

Mgt = II E '^+(^-'" ■ '^bT^^iW^t). 



N 



MT = Y.{0+\\J2h+irab,EN) 

ayib 



N 



(Ca ■ Tafc) (ffe - fab) - g CTq ' CT b 



(50) 
(51) 

(52) 
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N a^b 



N 



ajib 



M^ = J2{Q+\\J2K(rab,EN) 
N a^b 



a ■ ^ab){'^b ■ tab) — g^a ' 



rKllO.^)' 



M'P = E(0/ II E'^VK6,i^iv)^ {K - at) ■ [Tab X r+ab]}TlTl\\Ot), 



N ajib 



'^Tab 

R 



^'r = Y.^^1\\Y.Ki''-b,EN)—Vab ■ (o-a X Dfo + D„ X o-6)r;T^ 1 10+), 



AT 



Mgt^ = ^(0+||^/ioa.(rafe,-EjvK •«r6r^T|||0+), 
Mf^ = 5^(0+115^ /ioc.(r„6,i;ivKT^||0+), 

a=jLb 



N 



Mfo = E<0/ 1 1 E /^ol^a;., E^)v . v+tIt\ 1 10+) , 

N a^b 

Mpo = E(0/ II E '-§K{rab,E^)<T+ ■ [f X f+]r^r^||0+), 

AT a#6 

MGT = E(0/llE^+(^a6,i?iv)-a„.a,T>i^||0+), 
AT 05^6 

iR " 
r 



M^' = J2(^j\\J2Kirab,EN)- 
a^b 



N 



(fa • tab){crb ■ Tab) - -CTa ■ (Tb 



R 



M'nr = E(0/ II E Kirab,EM)^rab ■ (T„ - T(,Kr|||0,+), 



JV a^b 



2r 

iR 



M'rC. = E(0/ II E h'+irab,EN) — {iab ■ <TaCb ' CaVab ■ <T(,Kr|||0+), 



JV 



ai^b 



JV 



^KT. =E<'^/llE'^'+(^-«"^A^)?i'a6 • X Tfe + Ta X (T6)r>^||0+) . 

a#6 



In the above expressions, the neutrino potentials hi{rab, {En)) are defined as follows: 



h+{rab,{EN)) = 
hoirab, (En)) = 



R f dk 



+ 



1 



47r2 J u> \uj + Ai UJ + A2 
1 /■ / 1 1 



27r2£i 



a; V <^ + ^1 UJ + A-2 



2H{r,A)+r-^Hir,A), 



houj{rab, {En)) =h+- ARho, h'^{rab, {En)) = h+ + ARho, 



hR{rab,{EN)) = - 



A 



2 f R 



-ARh+ 



with 



H{r,A) 



277^ J (J (J 



,ik-] 



+ A' 



Aj = Sj + {En) -Ei, i = 1,2; A= {En) - {Ei + Ef)/2, 



(53) 
(54) 

(55) 

(56) 

(57) 
(58) 
(59) 
(60) 

(61) 

(62) 

(63) 

(64) 

(65) 

(66) 

(67) 



(68) 
(69) 

(70) 



(71) 
(72) 
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where rat is the distance between the nucleons a and b, and (-Ejv) is the average energy of the intermediate 
nucleus N. 

To derive the expressions for A and B shown in Tables 2 and 3 we have used the formulas: 

S-P p(3) 

^1 \^Stp) ~ ^^GT— —Xf 



rA(c^-p\—n (^ii+p\ ' —oi\t c^+^^,sp 

Ci (4^?) = 0, Ci {4%l) f = MaTel%^,x%; (73) 



{elz^) = Mgt {^i + 2/i^:^) ixF 1), 
(ey;^) = ^GT [m(x^^ - 1) + 2fi^-i{xF + 1)] , 

Zln{e^ti)-MGT4tix'H; (74) 



7.(3) 
Tl£_i_ 

9 A ' 



^ 9A 

W^2^(4^)--2zMgt4«^Xgt, 

/ NT r'^^ /I \ 

^^"^ ('^«) ^ = -4^Mgt4«^ (^gXGT + SXtJ • (75) 



(3) 



For all other arguments these nucleon matrix elements have zero values, except for 

Zi = 0) = MGT^i{XF - 1). (76) 

iSP T) 

We have calculated the numerical values of the integrated kinematic factors Aqi, Aq^ ' , i?oi, and 

(SP T) 

Bq^ ' for all the five nuclei of current experimental interest. We shall use them in the results shown below 
in Table 6 for the angular coefficient K. However, as we will focus in this paper mainly on the 01^2/3 decay 
of the ^^Ge nucleus, we give the values of these factors for this nucleus in Table 5, where we have used 

Q = E,-Ef- 2me = 2.039 MeV , (77) 

taken from Ref. [35], and the scaling factor for the neutrino potentials is 

R^roA^/^, ro = l.lfm. (78) 

The values of Aqq and Bos are of the order of 10"'*'' yr^^. Hence these values are not given in Table 5 and 
the terms with Aqq and i?03 can be safely neglected. 
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Table 5: The integrated kinematic A- and B-factors [in 10 ^^yr ^ for the 0+ 0+ transition 

of the 01^2/3 decay of ^^Ge. 



Am 


6.69 


-^01 


5.45 


Ao2 


1.09x10 


]3q2 


8.95 


Aqs 


3.76 


_Sq3 




^04 


1.30 




1.21 


^05 


2.08x10^ 


Ba5 


7.27 


^6 


1.69x10^ 




— 


Ao7 


1.05x10^ 




7.72x10-* 


Ao8 


6.59x10^ 


Bo8 


4.97x10^ 


Ao9 


4.14x10^ 


Bo9 


3.00x10^ 


ASP 
^00 


2.55 






aSP 
^01 


3.77 




2.73 


aSP 
^02 


1.18x10"^ 


nSP 
Bo2 


7.20x10"^ 


ASP 
^03 


1.27x10-'^ 


nSP 
-°03 


3.71x10"'* 


^01 


6.03x10 


Boi 


4.36x10 


^02 


1.50x10^ 


Bq2 


1.40x10^ 


^03 


7.67x10^^ 


Bo3 


7.16x10'^ 



We recall that the analytic expressions associated with the coefficients Cy^^ given in this section and 
the values of Aoi from Table 5 confirm the results of Ref. 5 . The analytic expressions associated with the 
coefficients fyiji^, ef!]lp: h ^^'^ values of Aq^'^''^\ Bqi, Sof^'^"* from Table 5 transcend the earlier 
work. 



3 Analysis of the electron angular correlation 
3.1 Qualitative analysis 

If the effects of all the interactions beyond the SM extended by the i^ms, which we call the "nonstandard" 
effects, are zero (i.e., all = 0), then K — Bqi/Aqi. Its values are given in Table 6 for various decaying 
nuclei. We will concentrate on the case of '^^Ge nucleus in the following. In this case the correlation ([2T|) is 
proportional to 1 — 0.81 cos9. (Note that in the limit of me/{Ei — Ej) ^ we have a+ + /3+ = 7+ + <5+ 
and K = 1.) Tables 2 and 4 show that the presence of the "nonstandard" parameters Cy^^, Ej^J^ or e^^ 
does not change the value of K and therefore the form of the angular correlation. The presence of any other 
parameter does change this correlation. From the fact that there are no contributions due to P-wave and 
recoil effects to the scalar nonstandard terms in the closure approximation (see Appendix A) , it follows that 
the values of ^02^1 ^03^' -^02^ j ^''^d -6^3^ are small and there are two additional "nonstandard" parameters 
that do not change significantly the form of the angular correlation, namely, efJp. 



Table 6: The values of angular correlation coefficient K for various decaying nuclei for the SM extended by 

the vmS. 





^«Ge 


«^Se 


i'^UMo 




*^«Xe 


K 


0.81 


0.88 


0.88 


0.85 


0.84 



Using Table 1 and taking into account the fact that are suppressed in comparison with |e^| by the 
factor rrii/me (the chiral suppression), we find the coefficient K and the set {e} of nonzero e^s that change 
the 1 — 0.81 cos 6* form of the correlation for the SM plus i^ms, see Table 7 (the lower two entries). They 
correspond to the following extensions of the SM: vms plus RPV SUSY [10 , i^ms plus right-handed currents 
(RC) (connected with right-handed W bosons 5 or LQs T5]). Hence, the angular coefficient K can signal 
the presence of these NP interactions. 



13 



Table 7: The angular correlation coefficient K for various SM extensions for decays of ™Ge. 



SM extension 


{e} 


K 






0.81 


i^M+RPV SUSY 




-1<K <l 






-1<K <l 



We remark here that in our earlier analysis [19] we had neglected the P-wave and recoil effects, which is 
not a good assumption. Our current study shows that these effects give significant contribution to the terms 
with Cyl ^ and e^^ . Hence, they have to be included in any realistic analysis of the data, as and when it 
becomes available. Including them, not only the model called vm+ RC but also the model vm+ RPV can 
essentially change the angular coefficient K from being 0.81 in the decay of the ^^Ge nucleus. Left-right 
symmetric models belong to the class z/m+ RC and we have studied these models in detail in section 4, where 
the correlations among the parameters if, T1/2 and either mw^j or C are worked out for the case 7^ 0, 

cosipi = considered in section 3.2. 

Note that the decay half-life and angular correlation do not give any bounds on the parameters e^^ and 
because the according expressions for A and B do not depend on them. 

3.2 Quantitative analysis 

Let us now consider some particular cases for the parameter space. We will analyze only the terms with 
^v^A corresponding nuclear matrix elements have been workd out in the literature. We use various 

types of QRPA model for the ^^Ge nucleus [2Ql|2T] as a test case. 

Using the case of \{m)\ — 0, which gives conservative upper bounds on and \e^\, the decay half-life 
is expressed from Eq. (pij) as 

ri/2 = ln2/r = (IMgtP^)^' . (79) 

From Eq. (|79|) . using Tables 2, 5 and the values of the nuclear matrix elements reported in Refs. [20ll21j . 
we have the following expressions for the half- life [in yr] for various choices of the parameters and 
[eyljl^l, taking only one parameter at a time: 

ri/2 - 1.1(1.3) X 10i2|^^:^|-2, Ti/2 = 3.2(4.0) X 10i2|^^;^|-2, (80) 
Ti/2 - 4.0(21) X 10i2|^^:^|-2, Ti/2 = 4.5(6.8) x 10i2|^^;^|"^ (81) 

ri/2 = 3.7(27) X W^\e}^ti\-', T,/^ = 1.0(9.7) x 10"|e^+^|-2. (82) 

Eq. (|5D|) corresponds to using the pnQRPA model with particle-particle strength parameter (;pp=1.02(1.06) 
PT] and Eqs (|5T|) - ([5^ correspond to using the QRPA model without (with) the p-n pairing [501 (note that 
the definitions of the nuclear matrix elements x'p and xn. in Ref. [20] differ from x'p and x'r in Ref. [S] by 
the factors 1/2 and A/{'meR), respectively). Comparing the numerical results in these equations, we note 
that the dispersion in the half-lifes is less marked for the coefficient However, the half-lifes involving 

the coefficients and |ey+^| show a very strong nuclear matrix element dependence. For the QRPA 

model worked out in [20j , it is not clear to us if this is due to a numerical artifact or the treatment of the 
isoscalar neutron-proton pairing. An important, and related point, is how to fix correctly the particle-particle 
strength of the nuclear Hamiltonian. Fixing the particle-particle pairing parameter, and varying it as done 
in [5T], leads to rather stable values for the half-life of ^^Ge nucleus. Clearly, these issues remain to be 
further discussed and clarified. A detailed discussion of these nuclear models will take us far afield from the 
main point of our paper. The theoretical uncertainty in the nuclear matrix elements [2l[36] plays an essential 
role in the numerical analysis. However, as we show below, the nuclear-model dependence of the angular 
coefficient K is rather modest. 
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The fact that the dependence of K on the nuclear matrix elements is much weaker than the uncer- 
tainty in Ti/2 from this source is illustrated in Table 8 for QRPA models [20l [21] for the assumed val- 
ues of the parameters: Imv^aI — l^y+^l = 5 x 10""^, |ey^^| = 5 x IQ-^. It is clear from Table 8 that 
measuring K with 10% accuracy (or better) produces useful experimental data that could be sensitive 
to the new physics. We note that for the parameters fJ^y^^ the angular coefhcient does not depend ac- 
tually on the nuclear matrix elements as it is seen from Tables 2, 3 (for \fj,\ — 0) and Eqs. (28), (31): 

K^iXFT if Boi/ \{xF T if Aoi] = BqiMoi ~ 0.81. 



Table 8: T1/2 and K for the fixed values of the parameters ky^^l for decay of ^^Ge for the case of 

|(m)| = in QRPA without (with) p-n pairing ^ [pnQRPA with ^^^=1. 02(1.06) [H]]. 





ImI^i^II = 5 X 10-^ 


K;^| = 5xio-^ 


=5x10-9 


\e'yXi\ = 5 X 10-^ 


^1/2/(10^^ yr) 


1.6(8.4)[0.44(0.52)] 


1.8(2.7)[1.3(1.6)] 


1.5(11) 


4.0(39) 


K 


0.81(0.81)[0.81(0.81)] 


0.81(0.81)[0.81(0.81)] 


-0.73(-0.73) 


-0.79(-0.87) 



Using the numerical results given above, the current lower bound T1/2 > 1.6 x 10^''' yr for the ^^Ge 
nucleus [37] yields the upper bounds on the parameters |m1/^^aI ^^'^ \^v^a\ shown in Table 9. The bound on 
kyl^l is stronger than the others shown in this table due to the relatively large values of the recoil and P- 
wave matrix elements in this case. The bounds on |eyi[^^| given in Table 9 are comparable with the bounds 
lentil < 4 X 10~9, \eyl^\ < 6 x 10^^ given in Ref. [H]- 



Table 9: Upper bounds on l^iy^^l, leJ/iJ^Al fo'' decays of ''^Ge for the case of = in QRPA. 



Nuclear model 




Kill 


\'^v-a\ 


ly+A, 
\W+a\ 


pnQRPA with .gpp=l. 02(1.06) ^ 


2.6(2.9) X 10"'' 


4.5(5.0) X 10-' 






QRPA without (with) p-n pairing [20, 


5.0(11) X 10-' 


5.4(6.5) X 10-' 


4.8(13) X 10-9 


7.9(25) X 10-' 



To be definite, we use the QRPA model without p-n pairing [20] in the following. The bounds given in 
Table 9 could be used for deriving the bounds on the parameters of the particular models (see section [2Tl) . 
For example, using Eq. ([5]) we have the following conservative constraints on the couplings of the effective 
LQ-quark-lepton interactions: 

• Consider a more general case of |(m)| 7^ 0, cosf/'i = 0, where the index i depends on a, (3 (as above, 
we take only one nonzero at a time). Using Tables 2 and 4 we have 

KA = DM^ +AD,\iii\\ (84) 



and 



KA^Di\ii\^ +D,\ei\^. (85) 



Hence, using Eq. (|79p we obtain 



\fi\^ = (Ai - A2if)/Ti/2, 

|e^P = (-A3 + A4if)/Ti/2 = 4|^^p, (86) 
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with the coefficients 



Ai = 

A3 = 



Di 
\Mgt?X 



A2 
A4 



C^ 

|A/gtPA,' 
Ci 



(87) 



where = CiD 
Using Eqs 



we have for e^+^ ^ 



(7.9 + IQK) X |e{;+:tl' = (5-1 - 6.3if ) x 



and for ej^l^ 7^ 



lAil^ = (7.7 + 10/0 X \el+_i? = (1.9 - 2AK) x 107ri/2 



(89) 



with Ti/2 in years. Fig. 1 shows the correlation among I1/2, K {left) and the correlation among 

kv+Al' ^1/2: -f'^ {fight) for the choice of a nonzero ey+^- Fig. 2 shows the same for the parameter Cyl^. It 
is clear from Figs 1 and 2 that the closer is X to 1 for the fixed value of ri/2, the weaker is bounded 
and stronger is bounded |ey^^|. The correlations among |eyj]^^|, 71/21 K will be used in the next section in 
the analysis of left-right symmetric models. 

Note that if several are nonzero in the considered model than the respective interference terms should 
be taken into account. 

• To extract |^|, |e(^|, c; in the general case of |(to)| 7^ 0, c; 7^ we need to analyze the data on 

at least two decaying nuclei. This analysis will be presented for the five nuclei already discussed in a 
forthcoming paper [38] . 



4 Electron angular correlation in left-right symmetric models 

The experimental bounds on the are connected with the masses of new particles, their mixing angles, and 
other parameters specific to particular extensions of the SM O HI [H [TOl [121 113] • To illustrate the kind of 
correlations that the measurements of ri/2 and the angular correlation coefficient K in the Qv2p decay would 
imply, we work out the case of the left-right symmetric models [52]. In the model SU{2)l x SU{2)ji x U{1) 
the parameters ry and A (see Eq. ([2])) are expressed through the masses and of the left- and 

right-handed W bosons and their mixing angle ^ |S]: 

ry = -tanC, X ^ {mw^ /mwa)'^ , (90) 

under the condition 

< fnwji- (91) 

Eqs. ([2]) and ([6]) and the relation [5] 

Ve^ = K'i (92) 

of the SU{2)l X SU{2)r x U{1) model yield 

el+i=X^Ue^Ve^., e^t^ = e^Ve, . (93) 

To reduce the number of free parameters, we assume the equality of the form factors of the left- and right- 
handed hadronic currents: 

9v = gv- (94) 

The small masses of the observable I'S are likely described by the seesaw formula that in the simplest case 
gives 

m, ~TO|,/Mi^, MH>mD, (95) 
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with the Dirac mass scale m/j (for the charged leptons and the Hght quarks niD ^ 1 MeV) and the mass 
scale Mji of right vms (in the majority of theories Mfj > 1 TeV). In the left-right symmetric models these 
scales arise usually from the two scales of the vacuum expectation values of Higgs multiplets [52]. In the 
seesaw mechanism, the values of the mixing parameters Vei (for i numbering light mass states) have the 
same order of magnitude as mD/Mn. In our discussion we use two rather conservative values (compare with 
Eq. (HI)) 

e = 10-^ 5 X 10^^ (96) 

for the mixing parameter 

e=\UM. (97) 

We recall that here the summation index i runs only over the light neutrino mass eigenstates (the summation 
over the total mass spectrum including also heavy states gives strictly zero due to the orthogonality condition 

From Eqs. ([901), ([Ml), dMl), and ([971) we have 

1 /2 

mwR=mwL{^/\(^vXA\) ' C = -arctan(|e^+^| /e) . (98) 

Using Eq. (pij) we note the approximate equality of mw^ aiid the mass of the observed charged gauge boson 
Wi (mvi/i=80.4 GeV [I]). 

The correlation among (C)i and T1/2 for the case of |(m)| 7^ 0, cosi/'i = (see section 3.2) 

is shown in Fig. 3 (4) for the two chosen values of e. The numerical results for these figures have been 
obtained using Eqs. (ISS)) and ([59]) . It is clear from Fig. 3 (4) that the closer is if to 1 for the fixed value 
of Ti/2 the stronger is the lower bound on mwn (the upper bound on Q). However this bound is weaker 
than the one mwa > 715 GeV, obtained from the electroweak fits [1]. There is still a more stringent bound 
> 1-2 TeV, obtained in Ref. [39] for the 0z^2/3 decay mediated by heavy Majorana neutrinos using 
arguments based on the vacuum stability [6] and additional theory input. We assume mwa > 1 TeV in the 
next figure. 

While experiments in the 0^2^ decay would measure the product of the quantities called A and the 
neutrino mixing matrix elements UeiVei in Eq. I|93p . collider experiments at the Tevatron and the LHC can, 
in principle, measure A by determining mwn- Assuming these logically independent possibilities, we plot 
the differential width ([5T|l vs. cos 6* in Fig. 5 for a set of values of |(to)| and mwn, taking e^^^ at a time 
and assuming e = 10^^. In this figure, we consider the values of |(m)|, starting from |(m)| < 0.03 eV up to 
I (m) I ~ 5 meV, covering two of three scenarios of neutrino mass hierarchies and mixing angles: normal and 
inverted mass hierarchies (see Ref. [ID] for a recent discussion and update). It is seen that the sensitivity of 
the electron angular correlation to the right-handed W^-boson mass tovFr increases with decreasing values 
of the effective Majorana neutrino mass |(m)|, as can be seen from Fig. 5 (right), where this correlation is 
shown for |(to)|=5 meV, 10 meV. 

In conclusion, we have presented a detailed study of the electron angular correlation for the long range 
mechanism of Oi/2/3 decays in a general theoretical context. This information, together with the ability of 
observing these decays in several nuclei, would help greatly in identifying the dominant mechanism underlying 
these decays. At present, no experiment is geared to measuring the angular correlation in Oi^2/3 decays, as 
the main experimental thrust is on establishing a non-zero signal unambiguously in the first place. We note 
that the running experiment NEMOS has already measured the electron angular distribution for the two 
neutrino double beta decay, and is capable of measuring this correlation in the future for the Oz/2/3 decay 
as well, assuming that the experimental sensitivity is sufficiently good to establish this decay [H]. The 
proposed experimental facilities that can measure the electron angular correlation in the 0i>2P decays are 
SuperNEMO [42], MOON ^43j, and EXO 44J. We have argued in this paper that there is a strong case in 
building at least one of them. 
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A 0u2P decay rate for scalar nonstandard terms 

The nucleon currents in the impulse approximation in the nonrelativistic form are used in this paper [321 134j . 
Keeping all terms up to order p/mp in the nonrelativistic expansion we have 

J+^p(x)^^r^5(x-r,)(Ff B,), - (99) 

a ^ 
J(^+^(X) = J2 - I-a) [g^'^igvla - 9ACa) + 9^''\gA<Ta,n " 9V " 5^^^)] , (100) 



2mp q^ + 2mp \ gv J 2mp q^ + 



(101) 



where = — p'^ is the 4-momentum transferred from hadrons to leptons, = +p'^; p^ and p'^ are 
the initial and final 4-momenta of a nucleon; is proton mass and TOtt is pion mass. 

(3) (3) 

We neglect the dipole dependence of the form factors Fg , Fp , gv, gA, gM on the momentum transfer 
and omit the zero argument of the form factors. Note that gv{0) = 1- 

Consider the pure SP case assuming (m) = 0. In terms of the combinations of hadronic currents 

J^^ = (F\J+\N)(N\J^+\I), J^^ ^ {F\J^^+\N)(N\J^\I), (102) 



i£+ = c/e. , (104) 



and the combinations 



r „ sf;'^(2y, Ix) s^-^ ly,2x) 

= ^ — ^ - ^ — ^, (105) 

sf, (2y, Ix) sf,(ly,2x) 
- ^ ' - ^ ' (106) 

of electron currents 

s^'^(2y, Ix) = e2(y)7A.(l T75)eJ(x), si^(2y, Ix) = e-2(y)7A(l - l5h^.el{x), (107) 
ei(x) = ep.s;(x), the matrix element is expressed as 



/2! V \/2 y (27r)3 2u; 

X ^ [m. {r^l^ - Jl_ll) + fc^ [Jl^l\^ - Jl^^l;)\ , (108) 

N 

where r = y — x. By using the identities 

s;;-«(ly, 2x) = s;:-^(2x, ly), si^(ly, 2x) = -s^^^{2^, ly), (109) 

the algebraic formula 

2{am ± bn) = (a + h){m ± n) + (a - h){m =F n), (110) 

the constant 



18 



and the neutrino potentials 



(27r)3 ijj UJ + Aj 



the matrix element (jlOSp is expressed as 



(112) 



(113) 



Each part of this matrix element is expressed as a sum of nonvanishing (indexed by n) and vanishing 
(indexed by c) terms, in the closure approximation: 



{M^p}n = I / d^dyT^iH, + H2) 

{M^p}, = I y dxrfyTjv(i/i - i?2) 

(Ai + + (A^ + i^pjF^_ + B.rFI + {B\ + SJpJf^i 



{M|p} 
{Af|p} 



2mp 



fci-^L) -(A2 + A2fi)i?^ + (4^- + i^^p)i?t + (Bi + syi?_ }, 



2mp 



d^dyTN { [H^i + i/. 



u)2 



-(hIi+hL) -{A2+A2R)Ei + {A'^' 



-{A\ + A\j,)E'_+B2rE+ 
A'J:^)E'1 + {B[ + B[^)E+] }, 



with 



The electron currents are defined as: 

P+ = \ ["(yx) ± u(xy)] , 
F^ = iK(yx)±u^(xy)], 

= iK''(yx)±w^nxy)], 



^5± = 2 K(yx) ± U5(xy)] 



^5± 



E± = F± + F. 



with 



5±, 



u(yx) = e2(y)e5(x), 
■u/^(yx) = e2(y)7''e!j;(x), 



K(yx) ±<(xy)] , 



[wr(yx)±u^''(xy)] 



■ 5±' 



U5(yx) = e2(y)75ei(x), 
U5(yx) = e2(y)757^ej(x), 



(yx) = -*e2(y)^^''ej(x), ^^(yx) = -*e2(y)75^'^''e^ (x 



The nucleon operator matrix elements are defined as follows: 

a = a + ap, b = b + bp, 



(114) 



(115) 



(116) 



(117) 

(118) 
(119) 



(120) 



(121) 
(122) 
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= 2G^£s, 


Air 


= — G° £sG+ 


— Lryepn+, 






A2 


= 2Gye5, 


A2B. 


= -G^e5G+ 


+ Gye'pB+, 






Al- 


= G°Aesal, 




= -G^\epBl_ 


f - Gy£s£''|_ 


APi _ 


-G^esP;, 




= G5ie5(T^, 




= G\e'pB]^j^ 


- Gye's 


aPi _ 




Af 


= i£ilkA\, 




— i£ilkA\p, 









BlR 

B2R 

3R ~ 



-G^esC^+G^'yepB^, 



GyS'pB^, 



-G\epBl_ - G°yesDl 



B\r — G\e'pBl 



G°ys'sDl 



r>Pi 

Bih = -G^.e'^P: 



G\esPl, 



(124) 



with 

B±=Bah±IaBb Bl^^alBh±Baal ^ P^h ± hPb- (125) 

Under the exchange of running indices a and b (i.e. x y), nuclear operators A, electron currents 
and and neutrino potentials Hi and H^i are even, while B, F^, and Hki are odd. 
The constants are defined as: 



G\ 



9v_ 
9A 



^V-A.i 



V-A 
^V+A,i 



G° = G(e = 0), 



G, 



^V-A 

-V-A.'. 



.V-A 

-V+A,i^ 



^0 _ 9V jj ^0 _ rr 

9A 



(126) 



p(3) , , p(3) 



'^S - —— Us-P,i + '^S+P.i) ^ 



5A ^ ' ' ' 9 A 

p(3) , , p(3) 



P I ^S-P _ S-P 
-P,i ^S+P.i 



9A ^ ^ 9A ^ ' 

Note that in the notations of Ref. [5] : 

< = u + M5, t'==w°'+u°'. (128) 

Since the nucleon recoil term behaves as an even parity operator while the neutrino momentum k 
and the recoil terms Ba, Ga, Da as odd ones, each of the Aj, k • A^, _Bj, k • Bj has a definite parity. The 
operators 

^1: ^3' ^4: ^3P' ^4Pi ^3' -°3P' 

r.B4p, r'A2fl, , (129) 



have even parity and the operators 



^IP' ^3P' ^4P,! ^IPi B2R, B\p, 

r • B4, r • Bf^, r'A2, r'4^ r'Afj^'^ , (130) 

have odd parity. The odd-parity operators do not contribute to the 0+ J+ transition in the case where 
both the electrons are in the 5'-wave state (the S ~ S case) with no de Broglie wave length correction (no 
FBWC). 

Using the definitions of neutrino potentials 

1 P 
h+ — —{Hi + H2), ho — — {Hi — H2), ft-ow = — {H^i — H^2), 
^ £21 £21 

= -^"4 (^^1 + ^^2) , Kf' = -^-^ (Hli - Hi,) , (131) 
^ £21 
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in the S ~ S case with no FWBC, Eqs. PT5)) . pTf)) are reduced to 

{MFp}n,s-s - / d^dyTr,h+ [A,F?^+ + (A^ + AP'i,)Fl^] , (132) 
{A^5p}c,5-5 = ^ / d^dyT^hoiBl + B[^)Fl, (133) 



+ 1 £Hi /■ dxdyT^-/.^f' (-^2Ki?i + A^inEl) , (134) 
{^^sp},,5-s = / d^dyT^Y^v ■ B4fl£;+, (135) 



where F are taken for x=0, y=0. 
For the 0+ 0+ transition we have 



Y.^X^^sp}s-s = 9\CfFL^ (136) 



N 



EE{^^|p}5-5 = 9l-^{C!^}cE^. (137) 

I N 



with 

Ci^ = (^/i+Ai), = (^/i'+f . B«), (138) 



where f = r/r and (X) = E E(0|ll^l|0|), with h = h{r,EN). 

i N 

In the S — Pi/2 case with no FBWC for the 0+ — > 0+ transition we have 



{MFpKs-p,,. - / d^dyTMh+ {Al,i,F^+ + Blj,Fl) , (139) 

Wsp}c^S-P,,, = ^ / '^^^y^^^o (^3K^5- + S^n) ' (140) 

i / d^dyT^ — h'^f' [-A2EI + {Al' + ^fA')^^+] , (141) 
{^^^5p}.,s-P,,, = J d^dyTr,ho^A\^El 

J dxdyT^^^/^;f' [-A^eL + (Ai^ + A^j,')E'l] . (142) 



l£2 

'2 m, 



2 



UleR 



The squared modulus of the matrix element (I113P , summed over the polarizations Sj of the electrons and 
multiplied by the phase space element ([55]) . yields the differential decay rate for the 0+ 0+ transition 

dT^Y. l^oTP£|df^o. = [A^ - Pi • P2Br] dno., (143) 

with aoiy being defined in Eq. ([36]). Here the coefficients are 

4 

A^^ = Y,\M,\^ (144) 

1=1 

= Re(MiM* + M*M2 + M3M4* + M*M4), (145) 
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with 



.2 



£2iR , X r 1 



m„R 


( < 


3 


\771e-R 


c 









2) C^R + '-^{Cinh^ ^ 
2) {{Cih - {C^}c - {C^r}c - {CtR}.) 



2meR 



{{C^}c + {CiR}c — 3{C^j?}c 



(146) 



M2 = at 



rrieR 

5 '^5R 



nipR 



nipR 



r 

2R 



c 



■(aZ)2 
M3 = at_i 
6 mei? 

Ml! 

2meR 



2R 6 Vwei? 

({C^}c + {CiR}c — '^{C^Rp} 

'2 



2 {{C^h - {C^h - {C^R^c - {<}c) 



(147) 



irieR 



5fl 



({C^}c + {C^r}c — 3{C^^}c 



(148) 



M4 = alii 



rrieR 



{^4r}c 



6 V ' 



C.^ - C, 



5R 



1 c 



6 TTleR 

[aZf 



({C^}c + {C^R^C — 3{C^J?}c 



2mei? 

where = yli(£2)^j (si) and the nuclcon matrix elements are 



(149) 



'-'3R — 


( 


■ Bsk) 




WLe r 




^3R — 


,mi i 


r+ • A 


3fl), {C3r}c = { — -hor- A3R 

nie r 


'^4R — 












A2), 





{Carf} 



'2R 
,iRrl 



^ h'r ■ BiR) 
2r 2E2 + 



(150) 



with r+ = y + X = 2i?f +. 

The terms in the first brackets in Eqs. (|146p - (|149p come from the S — S case, the terms in the second 
brackets come from the S — P1/2 case and in the third brackets there are the most important terms due to 
the Pi/2 - Pi/2 case and FBWC. 

Assuming now (m) 7^ for the dominant terms we have 

2 



Ml = all. 



nieR 



{C4:R}c 



'^^"^ {ci - Ct) ^ + i 



6me 



2R 6 \ meR 



c 



-2 {{Cth-iCth 



(151) 
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Ma = all 
£21^ , 



nipR 



M: 



6me 

3 = "l-l 
£2lR I £21 



c 



X 



2R 6 V rrieR 



+ 2 ({C2^},-{C5^}, 



TTLpR 



nie 



6 

£21^ / £21 



21(^2^-^5^) 



1 C 



2R 6meR 



X 



m„R 



6 



1 C 



6 mpi? 



({C^2^}c-{C5^}c) 



({C^2^}c-{C5^}c) 



(152) 



(153) 



(154) 



In the expressions for Adi, ...,M4, the terms with C are due to the inclusion of the P-wave in the elec- 
tron wave function and those with are from the inclusion of the nucleon recoil effect. In the closure 
approximation there are no contributions due to the P-wave and the recoil effects. Note that some of the 
subdominant terms should be taken into account in case of large cancellation among the dominant terms. 

B decay rate for vector nonstandard terms 

In this appendix we in general follow the derivation of Ref. ^ . However in addition to Ref. [5] we keep in 
our calculations the terms associated with the parameters Cy^^ and the pseudoscalar form factor. 

The nucleon currents in the impulse approximation up to order p/rup in the nonrelativistic expansion are 



(155) 



with Ca, D"^, given in Eq. PUT|) . 

In terms of Sl^u, Va^v, J^a (o^?/^ — L.R) 5 the matrix element 



i N 



(2^)3 uj 



may be expressed as 



The analogues of the Eqs. (C.2.11), (C.2.23), and (C.2.24) from Ref. [5j are as follows: 

{MH^Ac = {Mmjc = I y d^dyTN{Hi - H2) [(Xi + Xir)E^ + [Y^ + Yl„)E\_ 

{M^^}„ = {Mv+A{a)}^ = — /dxdyTjv { (-ff^i - H^2) 

me J 

(X3 + X,j,)Fl + Y^nF^_ + {Xl + X[R)Fi + {Yi + ^6';^)^^] + {Hii + Hi^) 
{Xi + Xij,)F^_ + {Yl + yIj,)FI+ + (Xf + Xl\)F^_ + (Fe'^- + YtR)FlA } , 



{M^^}, = {Mv+A{a)}c 



R 



j dxdyTN { (H^i + H^2) 



(156) 



RV = Co. E E (t^T^va + <a) , M^x' = {M^aU + {M^A%- (157) 



(158) 
(159) 



(160) 
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\xi + Xij,)Fl + [Yl + Yl^)Fl_ + (Xf + Xl\)Fl + {Yt + ^'Ij^s-] } , (161) 

with X = X + X^, Y = Y + Y^ . The operators X and Y arc defined in fS], except for the operator 
Y^j^ = —Y^j^ which is defined to remove the minus sign from the Eqs. (|160p and (|16ip : Xi = Xis, Yi = Yis- 
The additional operators are 

Fs^' = zG^e^e,.,P:^+ - G+Pi , Fg^' = -iGaSaSu^P^^ -G+pL, (162) 



with 



P:^+^<ylPl+K'yl (163) 



Under the exchange of running indices a and 6, nuclear operators X , electron currents Ej^ and F^ and 
neutrino potentials Hi and iJi^i are even, while Y , E^, F^, and Hki are odd. 
New constants are defined as: 



r,,-i}Lff^+A I \ p^-^V+A _ V+A 

9 A ^ ' 

The operators 

Y vP . vi \rPi. 

have even parity and the operators 



(164) 



Xl,X[j^, X[^, r-XsR, r^Xg'' 
YlYll vY.RyYil, (165) 



^iR\ Yir\^5R, r-Xs, r-Xfj^, r^Xll', r'X^J:^'; 

Y,R, Fg'^, r . Y3, r . Y^^^, r'Y^\ r%P^^ , (166) 

have odd parity. 

Using the definitions of the neutrino potentials from Eq. (|131l) and 

h^^ ^{H^i+H^2) (167) 
in the S ~ S case with no FBWC we have 

{M{^A}n,s-s ^ J d^dyTNh+{Xi+X[j,)E+, (168) 
{M^a}c.s~s = ^ / d^dyT^hoiY^ + Y,Pj^)EI, (169) 
{M^A}n^S-S = — / d^dyT^ho^ [(^3 + Xf«)P^ + (X^ + XS)P|] 
+ / d^dyTN^-^h'^f^ [Yij^F^^ + Yll,Fi^] , (170) 
{M^aUs-s d^dyT^hUYl + Y.^j^F^^ 

+ f d^dyTM—Kf' {Xi^Fl + Xi^^Fl) , (171) 
rrie J r 
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where E and F arc taken for x = y=0. 
For the 0+ 0+ transition we have 



N 



EE{<^} 



s-s - 9a 



i N 



Tile 



with 



■»'/? 



In the S — Pi 12 case with no FBWC for the 0+ — > 0+ transition we have 

{M{?4}„,s-p,/. = j d^dyTMh+Yi^El, 

e2iR 



{Mva}c,S-Pi,2 - 



d^dyTNhoYl^E\, 
{M^^}„,s-Pv2 = — / d^dyTr,hoo.{Xi^Fl + Yl^Fi_) 

'Trie J 

J d^dyTj-^h\f' [{Xf + X^k^)F^_ + {Yl^ + Yi'i^)F^^] , 



nieR . 



{M^^},,s-P,/, = ^J d^dyTNK{Xi^Fl+Yl^Fi+) 
+ £H1 [ d^dyrJ—Kr' [{Xf + XP}t)F'l + (Ff + Yri'')Ft 



The decay rate for the 0+ 0+ transition takes the form 



where the coefficients are 



i=l 



with 



iVi = ali_i 



Zi + -Z^iK 



XP 



m.-R 



^iR 



nipr 



nieR 



2 Z 



Y 
IR 



£21-^ 

2me 



{Zir}c 



c 



3 \meR 
{aZf 



-'^]{zl + zl^ + {zl^U- + t 



1 £91 R 



+ 



nieR 



N2 = all 



i{Z^}c + iZ^RF) 



z^ + z^p + —-z. 



4fl 



2R 3 me 



nier 



Zqr - 2^{Zj}c + {ZJr}c) 



6 



nieR 
2 



2 Z- 



Y 

IR ■ 



£2lR 

2m, p 
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TUeR 
6 TOp 



zi 



X 



■^Z^jip) 



„XP ^21 I ryX 



Z5R + {ZsrIc) 



yX 



{{Zt}c + {Z^R'^}c) 



r 

2R 



a_ 



Zt + 



£21 



Z^ 



ryXP 



3 VTLp 



r 

2R 



-li- 

3 \me 

{^3fi}c) 

-U- 

3 V Top 



(zf 



ryXP 



r 

6^ 



CZIr 



'2 HZ 



Zqr — "^Z^jj) 



(183) 



2me)R^Z^j, 



(184) 



:e2i(£2i -2me)i?^Zj;^ 



(185) 



where the terms in the first brackets in Eqs. p82p - (|185p come from the S — S case and the terms in the 
second ones come from the S — P1/2 case. The terms in the third brackets in Eqs. (|182P " P83p are the most 
important terms of those that come from the P1/2 — -P1/2 case and from the S — S case due to FBWC. The 
nuclear matrix elements are 



^IR 



Zl 



7^ 



{zYnh-i^^hor^-Y.n), 



vie 2R 



'2R 



-hn„,r ■ Y, 



6R/ 



Z, 

{zn 



(h',fVY^^), {ZUe - {^Kf^r^Yl^), 



■X 
4R 

■X 



,1 



{zt}e^{-Kerix, 

r 



{z^i}e^^-h:^%xl^^). 



,iR ri 



• 2r 2i?2 



h'v ■ Y, 



bR) 



The dominant terms give 

Ni = a% 



TUeR 



:Z 



AR. 



UleR 



2 U« 



2R 



Zi 



X 



nxeR 



-Z 



AR 



+ 


[-1 








\meR J ^ 2R 



1-1 



11 





- —z^ 


+ 




me 




zf 


+ —z^ 


+ 




nie 





1 fe 



-21 



-- — +2 Z 



tX 



tX 



U'-^-2]Z 



(186) 

(187) 
(188) 
(189) 
(190) 



that agrees with the Eq. (C.3.7) of Rcf. [S] taking into account the correspondence with their notations: 



Zl — Zl, 
ZJr = ZiR, 



Z3 - Z3, 

7X 



ZiR — Z5B, 



Zq — Zg, 
Z4 = Z5, 



(191) 



and the fact that Z2 is absent, as we have calculated only the leading contribution of the parameters e^. 
Recall that in Ref. [5] the pseudoscalar form factor is not taken into account. However the terms associated 
with this form factor do not contribute to the dominant terms (jl87p - (|190p . Note that in the expressions for 
A^i and A^2 given above, the terms with ( are due to the inclusion of the P-wave in the electron wave function 
and the ones with Z|^ are due to the nucleon recoil effect. We remark that some of the subdominant terms, 
like those with Z^, {Zf }c, {Z^^j^, {Zf and 

Z5RF' should be taken into account in the case of large 
cancellation among the dominant terms. The same is valid for the contribution due to the pseudoscalar form 
factor gAPa which yields corrections at about 10 % to the dominant terms. 
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C decay rate for tensor nonstandard terms 

The nucleon currents in the impulse approximation up to order p/vrip in the nonrelativistic expansion are 
used [31 [31], Jf^^A from Eq. ITODj) and 



0(,X - Taj ((,5'^ 5 - g )la+g g SkmnC^ak 



"F -e^"'"^ [(ffpfe5cr0 ~ gpOgakjTak + gprgcjs£rsk^ak\ 



rpk 



(3) 



2T, 



(192) 
(193) 



where, as before, = — is the 4-momentum transferred from hadrons to leptons, — + p'^, 
p^ and p'^ are the initial and final 4-momenta of a nucleon. We neglect the dipole dependence of the form 

(3) (3) 

factors and on the momentum transfer and omit the zero argument of the form factors. 
Consider the pure J^.i? case assuming (m) = 0. In terms of the hadronic currents 



J' 



{F\J2+\N){N\J^-^l\I), J^^^l^ = (F|4^^^|iV)(iV|jn/), 



I 



7^+ — JT ■ 7^+ 
'J L — '-^ei'Jv — A 1 



(194) 
(195) 
(196) 



and the leptonic tensors 



^a;..(2y,ix) t^p.(iy,2x) 

UJ + Ai UJ + A2 ' 



uj + Ai 



UJ - 



A. 



t^,^(2y,lx) <2^„(ly,2x) 
0J + Ai 

CAJ2y,lx) i2 (iy^2x) 



(197) 
(198) 
(199) 
(200) 



with the electron currents defined as 



tipuC^y, Ix) = e2(y)7a(l - 75)crp^e![(x), 
^aAM.^(2y, Ix) = e2(y)7„(l - 75)7AC^M->eJ (x), 
^^,ya(2y, Ix) = e2(y)o-p^(l - 75)7„e!j;(x), 
^LAa(2y, Ix) = e2(y)CT^^7A(l - 75)7„eJ(x) 



(201) 



the matrix element is expressed as 



m — 



N 



(202) 



For the electron currents we have the identities 

ti^,(ly,2x) 



-Ca(2y,ix) 



iiA;..(ly:2x) = i2^;,„(2y,lx) 



(203) 
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Using Eqs (|110p . (jllip . and (|112p . the matrix element (|202p is expressed as 



-^0^1/ — C'oi/ ^ ^ 



TO,; 



-M 



TOfi 



i N 

= {Mj, }„ + {My }c 



with nonvanishing (n) and vanishing (c) in the closure approximation parts: 
{M^}„ = R JdxdyTNiHi + H2) 

X "(C/i + + (U3 + U^r)FU + VifiF" + {V^ + V^r)F1 

{M^}^ = Rj d^dyT^iHi - H2) 



(204) 
(205) 



(206) 



(207) 



{M*;}„ = — / d^dyTN{H^i - H^2) V2bE- + {Ul + UIj,)FI' + {U^' + UI'j,)F!' 



+ {Hl, + HI,) {VI + Vl^)E^ + {U2 + U2r)F'.' + (C/| + Ui^F^ + (t/^^- + U:,=^)F[ 



•Oj 



R 



dxdyTjv(F^i + H^2) V2rE+ + {Ul + (71^)^0^ + {Ul' + C/^]^)F!; 



- HI2) [{VI + ^^4fl)^^~ + {U2 + ^2/^)^^ + {ui + ui^)FX + (t/;^ + c/;],)F^- 

+(;7^-''^ + c/8ii^)i^f 



Oj 



where the nucleon operators are 



U = U + UP, v = v + v^, 



Ul — ~2G°^{eTi + £T2)craCrb, U^jf — G"^{eTi + £T2)Pa+^ 

UiR = GUeT,+£T2)D%^iG\{eT,+eT2)T:+. 

U2 = 2iG\{e'T,+e'T,)cTaa,, C/^^^ = + 4J^"+' 

U2R = -tG°y{e'T^+e'TjD^:+ + G°Ae'T,+e'TjT^+, 

Ul - -G?,(eT,+£Tj<, Ui^^~iG\{eT,+eT2)er,kPit, 

3R = G'5i(£ti +eT2)C'^+ ~iGy{eTi + £T2)£yfc£'i- 

4i?, — *G'^(£ti + £t2)^<^+ ~ G'y(e^^ + e'rpJei-jkF'i'L 



c/1 



(7, 



f7, 



5fl. 



-iG\{e'rp^ + £7^2)^*^+ + Gv'^Ti^ijkDi'L 



(208) 



(209) 



(210) 



28 



Kr' = +iG°Ae'T,+e'^J{P^^+Pi\), 

t^8fl'' = -^G^(4,+4.K-^^V-G^^i0(4,^f++4.^^;), (211) 



with 



V^R = -G°yieTi+eT2)D'J_-iG°,{eTi+eT2)T^-, = -G^(eTi + et^P^'L, 

V^R = -GUe'Ti+eT2)D:_+iG°,{s'T^+e'r2)T:_, = + ^JP^'l, 

^^fl = -GV^ti + £T2)Q_ + «G"^(£T, + eT2)eijkDi''+ 

+ iG'yisT, + ST2)Ti. - iG^(eT, + eT2)eijknl, 

^3B. = (^Ti + £T2 )^i3kPiX , 

Vi = G°.(£^^+£^JaL-2zG^(£^^+4JKx<T,]\ 

Vlji = — G54(e^i + £^2)^*^- + «G?/(eyj + e^Jeijfe-D^*i_ 

F« = iG2i(4^+4j£i,fePi^, (212) 
Ti=T:h±Ian, rj^ = alTl±Ty„rj^^cTlXl±Xial, X = D,T,P. (213) 



Under the exchange of indices a and b, nuclear operators U, electron currents P+ and neutrino potentials 
Hi and H^i are even, while V, F- , and Hfej are odd. 
The new constants are defined as: 

7^(3) ^(3) 



The even parity operators arc 

C/i, C/f«, FC/2«, [/^, [/f;^,, [/i, f/f^, fc^C/^^, C/^^ ^/6^^^ fc*C/^i^^ 

V'aH, ^3, 1^3fl, k-V4K; (215) 

and the odd parity operators are 

UiR, fc'C/2) k''U2R, U^R, Ulji, k^Ui, k^U^^, U^r, k^Uj\ k^U^^ , 

k'ui^\ k'U^^^^; Fifl, F2ii, F^^, k.V4, k-Vf^. (216) 
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Using the definitions of the neutrino potentials from Eqs. (|13ip and (|167p . in the S — S case with no 
FBWC we have 

{M^}n.s-s = 2 / d^dyTNh+ [([/i + U[„)F°+ + (U^ + , (217) 

{^T }c,s-s = ^2iR J d^dyT^ho [V.^nFl + (V^ + Vfjl)Fl] , (218) 



rrie 

+ ^ / d^dyTN'-^h'_^r-YiRE+, 



'6R )^ + 



dxdyTN V^j^E+ 



+ - 



£21 



iR 

dxdyT/v — hg-r^ 



U2rF'_ 



•Oi 



TT3 F'J 



Tjij pOj , jjtjk rpjk 



where E and F are taken for x = y = 0. 
For the 0+ 0+ transition we have 



N 



EEW} 



s-s 



i N 



m„R 



{2Wrj, + {W,V},)E 



+ 1 



with 



1 1 tig / / tg / / tg 



In the S — Pi 12 case with no FBWC for the 0+ 0+ transition we have 
{M-}„,s_p,,, = 2 j dxdyTr,h+ (UIrF^^ + V^rFI) , 

{M™},,S-P,,, = £2ii? j dKdyTr^ho {UIrFI_ + V^rFX) , 

{M^Us-P,,, = — f dxdyTN Ho^UIrFI' 



/ dxdyT^'-^h'^f^ 



{M^}c,s-p,,, 



+ 



£21 



TTleR 



dxdyTN h^Ul^F'^ 



dxdyTN— h'^f' 



{U2 + U^r)FI^ + {U\' + U^}^)F'^ 



The decay rate for the 0+ 0+ transition takes the form 



[Al - Pi • paSj] d^o^, 



where the coefficients are 



2 



Bl = Re(OiO* + 0*02 + O3O: + O3O4) 



(219) 

(220) 

(221) 
(222) 



(223) 

(224) 
(225) 

(226) 

(227) 
(228) 

(229) 
(230) 
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with 



Oi = a*_i_-^ 



c 



TUpR 



2 W- 



u 

37?. 



nipri 

_2 



O2 - all 



nipR 



6 
2 

rrieR 



+ 



3(aZ)2 



nipR 



C 



6 



+ 


3(aZ)2 / 




mpR \ 



WXrf + 2{W^2Yf}< 



O3 = a\_^ 



vripR 



(W^4fi + W/}c) 



/ £21 



Top, 



£2\R I £21 



771 Q 



3 

3 TOe-K 2 



+ 


■ 3(aZ)2 / 




TOgi? \ 



W^4fiF + 7;{W2Vf} 



O4 = "-11 



e2i-R{W^iV^} 



TOgi? 

£2ii? / £21 

3 



+ 


■ 3(aZ)2 / 







(231) 



(232) 



,(233) 



,(234) 



where the terms in the first brackets in Eqs. (|23ip - (|234p come from the S — S case and the terms in the 
second ones come from the S — P1/2 case. The terms in the third brackets in Eqs. (|23ip - (|232p are the most 
important terms of those that come from the S — S case due to FBWC. Note that in the S — S case there 
is the contribution to Eqs. (|23ip and (|232p from the {Hi^i + H^2) combination in Eq. (|209p . Therefore the 
contribution from the Pi/ 2 — P1/2 case should not be taken into account. 
The nuclear matrix elements are 



'2R 



mi 



(235) 



Assuming now (to) 7^ for the dominant terms we have 



01 = 

02 = a*ii 



TOp/l 



Zf + 2< + -—{Wrn + {</}c) 

TTleR 



(236) 
(237) 
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£2l-R / £21 



rUe 



+ 

£2lR ( £21 



rrieR 

+ 2)(W^f^-{<}c-{<}c) 
2 



X 



+ 



me. 



2 «i^-{<}c-{<}c) 



4 C 



3 nieR 



4 C 



3 mpR 



+ 



+ 



(238) 



(239) 



Again, in the above expressions, the terms with C, are due to the inclusion of the P-wave in the electron 
wave function and the ones with and {X = U, V) are due to the nucleon recoil effect. 
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Figure 1: Correlation between the neutrino effective mass |(m)| (left) [ |ey^^| [right], the angular correlation 
coefficient K, and the half-life T1/2 for the Oz/2/3 decay of ^^Ge for the case cos -01 = 0. 




Figure 2: Correlation between the neutrino effective mass | (to) | (left) [|ey^^| (ri^/ii)], the angular correlation 
coefficient K, and the half-life T1/2 for the Oj^2/3 decay of ^^Ge for the case cos'i/'i = 0. 
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Figure 3: Correlation between the right-handed M^-boson mass itlwri the angular correlation coefficient 
K, and the half-life T1/2 for the 01/2(3 decay of ""^Ge for the case cos'f/'i = and e = 10^^ (Isft) and for 
e = 5 X 10"'^ (right). 




Figure 4: Correlation between the mixing parameter C, the angular correlation coefficient K, and the half-life 
Ti/2 for the Oz^2/3 decay of ™Ge for the case cos V-i = and e = 10"^ (left) and for e = 5 x lO""^ (right). 
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Figure 5: Left: Differential width in cos^ for the 01/2/3 decay of ^^Ge for a fixed value of e = 10~^ and 

|(m)| = 20,30 mcV. The straight and dotted Uncs correspond to m,WR = 1 TcV, oo, respectively (the 
latter is the conventional case of the light Majorana neutrino exchange mechanism). Right: The same as 
the left figure but for smaller values of |(m)| = 5,10 meV. In addition, the dashed lines correspond to 
mwR = 1-5 TeV. 
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